LOCAL EXPONENTIAL H 2 STABILIZATION OF A 2 x 2 
QUASILINEAR HYPERBOLIC SYSTEM USING BACKSTEPPING * 
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Abstract. In this work, we consider the problem of boundary stabilization for a quasilinear 
2x2 system of first-order hyperbolic PDEs. We design a new full-state feedback control law, with 
actuation on only one end of the domain, which achieves H 2 exponential stability of the closed- 
loop system. Our proof uses a backstepping transformation to find new variables for which a strict 
Lyapunov function can be constructed. The kernels of the transformation are found to verify a 
6JQ[ Goursat-type 4x4 system of first-order hyperbolic PDEs, whose well-posedness is shown using 

the method of characteristics and successive approximations. Once the kernels are computed, the 
stabilizing feedback law can be explicitly constructed from them. 
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1. Introduction. In this paper we are concerned with the problem of boundary 
stabilization for a 2 x 2 system of first-order hyperbolic quasilinear PDEs, with ac- 
tuation at only one of the boundaries. The quasilinear case is of interest since many 
relevant physical systems are described by 2 x 2 systems of first-order hyperbolic 
quasilinear PDEs, such as open channels [TB I fT7 l I18|. transmission lines [6], gas 
flow pipelines [14] or road traffic models [12| . 

This problem has been considered in the past for 2x2 systems [T3| and even 
n x n systems [24] . using the explicit evolution of the Riemann invariants along the 
characteristics. More recently, an approach using control Lyapunov functions has 
been developed, for 2 x 2 systems [2] and n x n systems [3J. These results use only 
static output feedback (the output being the value of the state on the boundaries). 
However, they do not deal with the same class of systems considered in this work 
(which includes an extra term in the equations); with this term, it has been shown 
in [1] that there are examples (even for linear 2x2 system) for which there are no 
control Lyapunov functions of the 'diagonal" form J z(x,t) T Q(x)z(x,t)dx (see next 
section for notation) which would allow the computation of a static output feedback 
law to stabilize the system, even if feedback is allowed on both sides of the boundary. 

Several other authors have also studied this problem. For instance, the linear case 
has been analyzed in |38) (using a Lyapunov approach) and in |25| (using a spectral 
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approach). The nonlinear case has been considered by [8j and [15] using a Lyapunov 
approach, and in [22], [57], and |TT] using a Riemann invariants approach. 

The basis of our design is the backstepping method [20] ; initially developed for 
parabolic equations, it has been used for first-order hyperbolic equations [25], de- 
lay systems [2T] . second-order hyperbolic equations [ST], fluid flows [34], nonlinear 
PDEs [35] and even PDE adaptive designs [32]. The method allows us to design a 
full-state feedback law (with actuation on only one end of the domain) making the 
closed-loop system locally exponentially stable in the H 2 sense. The gains of the feed- 
back law are the solution of a 4 x 4 system of first-order hyperbolic linear PDEs, whose 
well-posedness is shown. The proof of stability is based on [3] ; we construct a strict 
Lyapunov function, locally equivalent to the H 2 norm, and written in coordinates 
defined by the (invertible) backstepping transformation. 

The paper is organized as follows. In Section [2] we formulate the problem. In 
Section[3]we consider the linear case and formulate a backstepping design that globally 
stabilizes the system in the L 2 sense. In Section 0] we present our main result, which 
shows that the linear design locally stabilizes the nonlinear system in the H 2 sense. 
The proof of this result is given in Section [5l We finish in Section [6] with some 
concluding remarks. We also include an appendix with the proof of well-posedness of 
the kernel equations, and some technical lemmas. 

2. Problem Statement. Consider the system 

z t + A(z,x)z x + f(z,x) = 0, x G [0, 1], t £ [0, +oo), (2.1) 

where z : [0,1] x [0,oo) — > M 2 , A : R 2 x [0,1] -> .M 2 , 2 (R), f : M 2 x [0,1] -> M 2 , 
with -M 2 ,2(R) denoting the set of 2 x 2 real matrices. We assume that A(z, x) is twice 
continuously differentiable with respect to z and x, and we assume that (possibly 
after an appropiate state transformation) A(0,x) is a diagonal matrix with nonzero 
eigenvalues Ai(x) and A 2 (a;) which are, respectively, positive and negative, i.e., 

A(0,aO=diag(Ai(x),A 2 (z)), Ai(x) > 0, A 2 (x) < 0, (2.2) 

where diag(Ai,A 2 ) denotes the diagonal matrix with Ai in the first position of the 
diagonal and A 2 in the second. 

We also assume that f(0,x) = 0, implying that there is an equilibrium at the 
origin, and that / is twice continuously differentiable with respect to z. Denote 



f2l(x) / 22 (x) 



(2.3) 



and assume that € C 1 ([0, 1]). 

Denoting z = \z\ z 2 ] T , we study classical solutions of the system under the fol- 
lowing boundary conditions 

Zl (0,t)=G (z 2 (0,t)), z 2 (l,t) = U(t), te[0,+oo) (2.4) 

which are consistent (see [28]) with the signs of (|2.2|) . at least for small values of z. 
We assume that Go (x) is twice differentiable and vanishes at the origin. In (|2.4[) , U (t) 
is the actuation variable, and our task is to find a feedback law for U (t) to make the 
origin of system (|2.ip . (l2.4p locally exponentially stable. 

Remark 1. The case with / = in (|2.ip was addressed in [2] and [3] by using 
control Lyapunov functions to design a static output feedback law; this approach has 
been shown to fail in [1] for some cases with / ^ 0, at least for a "diagonal" Lyapunov 
function of the form J Q z T (x,t)Q(x)z(x,t)dx. 
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3. Stabilization of 2 x 2 hyperbolic linear systems. Next, we present a 
new design, based on the backstepping method, to stabilize a 2 x 2 hyperbolic linear 
system; this procedure will be used later to locally stabilize system (|2.1[) . (|2.4p . 

Consider the system 

w t = Z(x)w x + C{x)w 1 x e [0, 1], t e [0, +oo), (3.1) 

where w : [0, 1] x [0, oo) -> R 2 , E, C : [0, 1] -> A< 2 , 2 (M), where the matrices £ and C 
are respectively diagonal and antidiagonal, as follows: 

=W-(- a W e£))' CW =(j*) t')' (3 ' 2 ' 

where ci(x),C2(x) areC([0,l]) and ei(x),e2(a;) areC 1 ([0,l]) functions, verifying that 
ei (x) , e% (x) > 0, and with boundary conditions 

u(0, i) = gru(0, i), u(l, i) = U(t), (3.3) 

where g € R and the components of u> are ti) = [u v] T . Our objective is to de- 
sign a full-state feedback control law for U(t) to ensure that the closed- loop system 
is globally asymptotically stable in the L 2 norm, which is defined as t)||^2 = 

\J Jo ( m2 (£; t) + w2 (C: *)) d£- There are two cases, depending on whether q in (I3.3[) is 
nonzero or q = 0. We first analyze the first case, thus assuming q =/= 0. 

3.1. Target system and backstepping transformation. Our approach to 
designing U(t), following the backstepping method, is to seek a mapping that trans- 
forms w into a target variable 7 with asymptotically stable dynamics as follows: 

7t = £(x) 7x , (3.4) 

with boundary conditions 

a(0,t) = q/3(0,t), 0(l,t)=O, (3.5) 

where the components of 7 are denoted as 

7(1, t) = [a(x,t) P{x,t)] T . (3.6) 

System (|3.4II . (|3.5p verifies the properties expressed in the following proposition. 

Proposition 3.1. Consider system {3.$ , i3. 5\) with initial condition 70 G 
L 2 ([0, 1]). Then, for every A > 0, there exists c > such that 

\h(;t)\\L> <ce- At || 7 o||L2. (3.7) 
In fact, the equilibrium 7 = is reached in finite time t = tp , where tp is given by 

"-jf(sb + ^)* <3 - 8) 



Proof. Define 



4 



J.-M. CORON, R. VAZQUEZ, M. KRSTIC AND G. BASTIN 



where A, B, n > will be computed later. Select: 

Vi= f j T (x,t)D(x)j(x,t)dx. (3.10) 



Notice that \fV\ defines a norm equivalent to ||tC"» ^)||i 2 - Computing the derivative 
of V\ and integrating by parts, we obtain 

Vi = - / j T (x, t) (D(x)12(x)) x j(x, t)dx + [ 7 T (x, t)D(x)V{xyi(x, t)] I , (3.11) 



> 0, (3.12) 



where we have used that S(a;) and D(x) commute. Since 

Ae-^ x 
Be^ x 



(D(x)E(x)) x =fx 

and, on the other hand, 

[j T (x, t)£>(z)£(*)7(z, *)] I = -Aa 2 (l, tje"" — (B — q 2 A)j3 2 (0, t), (3.13) 

choosing B = q 2 A + A 2 , A = A 2 e^, and fi = A x e, where e = max xe [ ,i] {77^7 > i^J }> 
we get that (D(a;)E(a;)) a , > Ai£>(a;), therefore: 

^ < -AiVi - A 2 (a 2 (M) + /3 2 (0,i)) , (3.14) 

where Ai,A 2 > can be chosen as large as desired. This shows exponential stability 
of the origin for the 7 system. 

To show finite-time convergence to the origin, one can find the explicit solution 
of (|3.4p as follows. Define first 

M*) = I -7Z7*> M*) = I ~7ts<%, (3-15) 



61 (0 ^ ' Jo *2® 

noting that they are monotonically increasing functions of x, and thus invertible. Note 
that the components of 7 verify the differential equations 

a t = -ei(x)a x , (3.16) 

Pt = e 2 (x)p x , (3.17) 

which can be rewritten as follows 

-a{cj>l x (x),t) + — a (^\x),t) = 0, (3.18) 



: ^ 2 - 1 (x),t) - /3(<j>-i(x),t) = 0. (3.19) 



The solution of these equations is a(x,t) = F a (<f>i(x) — t) and ft(x, t) = Fp(<f>2(x) + t), 
where F a and Fp are arbitrary functions. Now, if ao(x), Po{%) are the initial condition 
for the states, one obtains F a {x) = ao(</>j" 1 (a;)) (valid for < x < </>i(l)) and Fp(x) = 
j3o((f>2 (x)) (valid for < x < </> 2 (l)). Using the boundary conditions (|3.5[) one finds 
the remaining values of F a and Fp, and thus the solution of the system, as follows: 

a(M) = { a ° ^t^lV^ I > J SI' ( 3 ' 2 °) 
I qp{0,t - 4>i(x)) t>4>i(x), 

W'V-X t > 2 (i) - 4> 2 {x), [A - n) 
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Thus, after t = tp, where 

"=*' (i)+ « i) =r(^) + ^))«' (3 - 22 » 

one has that a = f3 = 0. □ 

3.2. Backstepping transformation and kernel equations. To map the orig- 
inal system (|3.ip into the target system (|3.4p , we use the following transformation: 



7 = w- K(x,£)w(£,t)d{, (3.23) 
Jo 

where 

{ ~ I K vu (x,0 K vv (x,0 J { > 



is a matrix of kernels. Defining 



0„=( J J ]. (3.2.-,) 



the original and target boundary conditions (respectively (|3.3p and (13. 5)) ) can be 
written compactly (omitting dependences in x and i) as 

u>(0, t) = Q Q w(Q, t), w(l, t)=(p\ 7(0, t) = Q o7 (0, t), 7 (1, t) = 0. (3.26) 

Introducing Q3.23P into (|3.4p , applying (|3.ip , integrating by parts and using the bound- 
ary conditions, we obtain that the original system (|3.ip is mapped into the target 
system (|3 .4[) if and only if one has the following three matrix equations: 

= C{x) + T,(x)K(x, x) - K(x, x)Y,(x), (3.27) 

= X(x)K x (x, + K^x, 0E(0 + K(x, £)£'(0 - 0<7(0, (3.28) 

= J ftT(x,0)S(0)g o . (3.29) 

Expanding Q3.27p . we get the following kernel equations: 

ei(x)K™ + e x (0Kp = -e'^K™ - c 2 (£)iT"\ (3.30) 

ex(x)iq? - e 2 (0*T = e' 2 (£)K™ - Cl (£)#"«, (3.31) 

e 2 {x)K™ - d(0^r = e i(0^" + c 2 (0K™, (3.32) 

£2(2;)^ + C2 (0*r = -4(0^™ + ci(e)^", (3.33) 
with boundary conditions obtained from (|3.28l) - (|3.29l) 

X"«(x,0) = ^- K ™( a; ,0), (3.34) 

K™{x,x)= f^fv (3.35) 

X™(z,z) = ^ (3.36) 

ei{x) + e 2 {x) 

K vv (x, 0) = ^j^-K va (x, 0). (3.37) 
e 2 (0) 
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The equations evolve in the triangular domain T = {(cc,£) : < £ < x < 1}. 
Notice that they can be written as two separate 2x2 hyperbolic systems, one for K uu 
and K uv and another for K vu and K vv . 

By Theorem I A. 1 1 ( see the Appendix), one finds that, for q ^ 0, under the assump- 
tion that ci(x), C2(x) are C([0, 1]), £i(x), £2(1) are C 1 ([0, 1]) and that ei(x), e 2 (x) > 0, 
there is a unique solution to (|3.30|l — (|3.37|1 . which is in C(T). 

3.3. The inverse transformation. To study the invertibility of transformation 
(|3.23p . we look for a transformation of the the target system f|3 .4[) into the original 
system (|3.ip as follows: 



w{x, t) = 7 (x, t) + / L{x, t)d£, (3.38) 
Jo 

where 

L(x €\ - ( La " M \ (3 39) 

Introducing (|3.38p into (|3.ip , applying (|3.4p , integrating by parts and using the bound- 
ary conditions, we obtain as before a set of kernel equations: 

e x (x)LT + ei(0^r = -t'i(Z) Laa + ci(x)L^ a , (3.40) 

e x {x)Lf - e 2 (0Lf = e' 2 ^)L^ + c x {x)L^ , (3.41) 

e 2 (x)L? a - £l (0if = e[(0L^ - c 2 (x)L aa , (3.42) 

e 2 (x)L^ + C2 (0Lf = -e' 2 (0L w - c 2 (x)L^, (3.43) 
with boundary conditions 

L™(x,0) = -^L°f>(x,0), (3.44) 

L^( I , ar ) = _^L_, (3.45) 



ei(x) + e 2 (x) 

L^(*,0) = ^L^,0). 
e 2 (0) 



(3.46) 
(3.47) 



Again by Theorem I A. II (see the Appendix), one finds that there is a unique 
solution to these equations, which is C(T). 

3.4. Control law and main result. From the transformation (|3.23p evaluated 
at x = 1, one gets 

U= [ /T»(l,0u(£, m+ I K™(l,0v(t,t)dt (3.48) 
Jo Jo 

With control law (|3.48[) the following result holds. 

Theorem 3.2. Consider system i3.1\) with boundary conditions i3.3\) , control 
law and initial condition Wq £ L 2 ([0,1]). Then, for every A > 0, there exists 

c > such that 

|«,t)|| L 2 < ce^H^oH^- (3.49) 
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In fact, the equilibrium w = is reached in finite time t = tp, where tp is given by 

Proof. Since the transformation Q3.23[) is invertible, when applying control law 
(|3.48p the dynamical behavior of (|3.ip is the same as the behavior of (I3.4[) . which is 
well-posed from standard results and whose explicit solution and stability properties 
we know from Proposition 13.11 Thus, we obtain the explicit solutions of w from the 
direct and inverse transformation, as follows: 



,{x, t) = 7* (as, t) + f L{x, 07* (6 
Jo 



(3.50) 



where 7*(x, i) is the explicit solution of the a, (3 system, given by (|3.20|) - (|3.21|) . with 
initial conditions: 

l0 (x)=w (x)- [ K(x,0wo(0d£. (3.51) 
Jo 

In particular, we know that 7 goes to zero in finite time t — tp, therefore w also 
shares that property. Finally, since the origin of the 7 system is L 2 exponentially 
stable with an arbitrary large exponential decay rate, we conclude, using the inverse 
transformation, that the origin of the w system is also L 2 exponentially stable with an 
arbitrary large exponential decay rate. Equation (|3.49p follows by using the inverse 
and direct transformations to relate the L 2 norms of w and 7 (using the fact that the 
kernels of the transformations are continuous, and thus bounded, functions). □ 

3.5. The case q = 0. If the coefficient q is zero in (|3.3[1 . the method presented 
in the paper is not valid since (|3. 341) would require the value of one of the control 
kernels to be infinity at the boundary of the domain T . Similarly, if the coefficient 
is close to zero one still gets very large values for the kernels close to the boundary, 
resulting in potentially large control laws. 

The method can be modified to accommodate zero or small values of q by setting 
a slightly different target system (|3.16|) - (|3.17D . as follows: 

at = -ei(x)a x + g(x)p(p,t), (3.52) 
Pt = e 2 (x)f3 X) (3.53) 

where g(x) is to be obtained from the method; regardless of the value of g(x), this is 
a cascade system which is still L 2 exponentially stable and converges in finite time by 
the same arguments of Proposition 13. 1[ since now, using the same Lyapunov function 
V\ defined in (|3.10[) . we obtain 

V x = - [ 7 T (x, t) {D{x)V{x)) x 7 (x, t)dx + [^(x, t)D(x)E(xMx, t)] J 
Jo 

+2/3(0, t) [ a{x,t)A^—g{x)dx, (3.54) 
Jo eiW 

The new term (which is the last one) can be controlled by slightly modifying the 
coefficients of D(x) in the proof of Proposition 13.11 obtaining the same result as 
before. 

The kernel equations resulting from the transformation are still the same (|3. 301) - 
with the same boundary conditions ([3T55 ) -([337 p for K uv ', K vu , and K vv 
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(which reduces to K vv (x,0) - 
boundary conditions for K uu : 



when q = 0), but one obtains an undetermined 
K uu (x,0) = h{x), (3.55) 



where h(x) can be chosen as desired; by choosing at least a continuous function, one 
can apply Theorem I A. 1 1 and thus the kernel equations are well-posed. After h(x) has 
been chosen and the kernels have been computed, one obtains the value of g(x) as 



g{x) = qei(0)h(x) - e 2 (0)K uv (x, 0). 



(3.56) 



Invertibility of the transformation follows as before, thus one obtains the same 
result of Theorem l3.2l The non- uniqueness in p.55[) gives the designer some freedom 
in shaping the input function g{x) from /3 to a. Also note that this has no impact in 
the feedback law as the kernels K vu and K vv (which are the ones appearing in (|3.48p ) 
are uniquely defined and independent of the non- unique K uu and K uv . 

4. Application of the linear backstepping controller to the nonlinear 
system. We wish to show that the linear controller Q3.48P designed using backstep- 
ping works locally for the nonlinear system, in terms that will be made precise. 

For that, we write our quasilinear system (|2.ip in a form equivalent (up to linear 
terms) to (|3.ip . Define 



ipi(x) = cxp 



x fn(s) 

A 1 (s) 



ds 



<p 2 (x) = exp - 



x Mg) 

A 2 ( S )' 



We obtain a new state variable w from z using the following transformation: 



w(x, t) 
so that 



u(x, t) 




v(x, t) 





ipi{x) 

ip 2 (x) 



zi(x,t) 
z 2 (x,t) 



= $(x)z(x,t), 



(4.1) 



(4.2) 



z(x, t) 



1 



tp 1 (x) 




(p 2 (x) 



w(x,t) = $ 1 (x)w(x,t). 



It follows that w verifies the following equation: 

w t + A(w, x)w x + f(w, x) = 0, 

where 

A(iu,x) = ^(x)A(^' 1 (x)w,x)^ 1 (x), 



f(w,x)=$(x)f($ 1 (x)w,x) + A(w,x) 



hi(x) 
A x {x) 







/22(S) 

A 2 {x) 



w. 



(4.3) 

(4.4) 

(4.5) 
(4.6) 



It is evident that A(0,a;) = $(x)A(0, x)^ 1 (x) = A(0,x) and that f(0,x) = 0. Also 

df(w, x) 



C(x) = - 



dw 



w=0 



-h 2 (x) 

-f2l(x) 



(4.7) 
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Thus, it is possible to write (|4.4p as a linear system with the same structure as p. II) 
plus nonlinear terms: 



where 



and 



w t - T,{x)w x - C(x)w + A NL (w, x)w x + /nl(w, x) = 0, (4.8) 



E(x) = -A(0,x), (4.9) 



A NL (w,x) = A(w,x) + 'E(x), f NL (w,x) = f(w,x) +C(x)w. (4.10) 



Computing the boundary conditions of 



formation (|4.2p . and defining q 



dGo(v) 



v=0 



by combining (|2.4p with the trans- 
and Gnl{v) — Gq(v) — qv, one obtains 



w(0, t) = qv(0, t) + G NL (v(0, t)), v(l,t) = U(t), 



(4.11) 



where U(t) — (p2{l)U(t). In what follows we will consider the case q ^ 0; the case 
q = is analogous (see Remark [5]) . 

Notice that the linear parts of (|4.8p and (|4.1ip are identical to (|3.ip and (|3.3p . 
and that the coefficients C(x) and S(x) verify the assumptions of Section [31 Also, it 
is clear that the nonlinear terms verify Ajvl(0, x) — 0, /atl(0,x) = d (^ L (0, x) = 0, 
andGAT L (0) = ^(0) =0 

Therefore, we consider using the feedback law: 



U 



K™(l,Z)u(Z,t)dt- 



(4.12) 



which implies, in terms of the original z variable: 



U 



^a(l) 



K vu (l,Z)z 1 (t,t)tp 1 (£)<%+ / A- W (l,0*2(e.*)^(0* ,(4.13) 



where the kernels are computed from ()3.30[) — f!3.37|) using the coefficients C(x) and 
S(x) from g3) and (gU). 

Next, we show that the control law (|4.13p . which is computed for the linear part of 
the system, asymptotically stabilizes the nonlinear system, although locally. However, 
the right space to prove stability of the closed- loop system is H 2 , instead of the space 
L 2 that was used in Section [3] for the linear system. 

Denoting: 



go 



,G{z) = G (z 2 ),qi 



,k(x) 



the boundary conditions of the closed loop system would be written as: 



(4.14) 



qZz(0,t) = G(z(0,t)), qjz{l,t)= / fc r (0z(£,t)d£ 



(4.15) 
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A necessary condition for system (|2.1[) with boundary conditions (|4. 15[) to be well- 
posed in the space H 2 is that the initial conditions verify the corresponding second- 
order compatibility condition. These are 

= G(z o (0)) - gfcoCQ), (4.16) 

o= / fc T (£M£K-^o(i), (4.17) 
Jo 

= G'(z o (0)) (A(z (0), 0)4(0) + /(z„(0), 0)) 

~ql (A(«b(0), 0)4(0) + /(z„(0), 0)) , (4.18) 

fc T (£) (A(z (0>0*o(£) + /(*o(0,OK 
-qf (A(z (l), 1)4(1) + /(z (l),l)). (4.19) 

While (|4.16p and (|4.18p are natural compatibility conditions, the conditions (|4.17l) and 
(|4.19p are artificial (since they show up due to the feedback law that has been designed) 
and rather stringent, as they require very specific values of the initial conditions. Thus, 
we modify our control law in a way that, without losing its stabilizing character, does 
not require any specific values in the initial values beyond the natural conditions 
(|4.16p and (|4.18p . The modification in the boundary conditions consists in adding a 
dynamic extension to the controller as follows: 



z(0, t) = G(z(0, t)), qi z(l, t) = / k 1 (£)*(£, t)d£ + a(t) + b(t), (4.20) 

Jo 

where a(t) is one of the states of the following system: 

a = -dia, b = -d 2 b, (4-21) 

where the constants d\ and d 2 can be chosen as desired with the only conditions that 
di, d 2 > and d\ ^ d 2 . It is evident that with positive values of these constants, (|4.2ip 
is always stable. The initial conditions of a(t) and b(t) are an additional degree of 
freedom that can be used to eliminate the compatibility conditions (|4.17p and (|4.19p . 
With the modification of the control law, these compatibility conditions are now 

l 



0= / fc J (0^o(?)^ + a(0) + &(0)-9i^(l), (4.22) 

k T (0 (A(3b(0. (0 + /(*>(£)> 0) d t ~ dMO) - d 2 b(0) 
-qf (A(z (l), 1)^.(1) + f(zo(l), 1)) • (4.23) 



Call 



l 



Pi{zo) = q{ Zo(l) - / k 1 (OzoiOd?, (4.24) 
Jo 

P 2 (z ) = ql (A(z (l), l)2d.(l) + /(*b(l), 1)) 

k T (0 (A(z (0,0z 0x (0 + f(zo(£), 0) (4.25) 



Selecting 



a(0) = -^M±^M ; 6(0) = ^M+^Cfo) f (4 26) 
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the compatibility conditions are automatically verified. 

We are now ready to state our main result. Define the norms ||z(-,i)|| ff i 
||«(-,i)|| L 2 + \\z x (-.t)\\ L 2 and IK,*)!!^ = ||«(-,t)|| HI + \\z xx (-,t)\\ L 2. 

Theorem 4.1. Consider system i2.1\) and with boundary conditions J^.l 

and initial conditions zq = [zq 1 zq 2 ] t € H 2 ([0, 1]), and a(0) and 6(0) verifying {4- 



with the kernels K vu and K vv obtained from f3.30\) - [3.37)) where the coefficients C(x) 
andT.(x) are computed from |^.7| ) and |^.ff[ ). Then, under the assumptions of smooth- 
ness for the coefficients stated in Section [H for every A > 0, there exist 5 > and 
c > such that such that, if \\zq\\jj2 < 6 and if the compatibility conditions \4- 16ty and 
{4 -Iffy are verified, then: 

M',t)\\h + «*(*) + b2 W < ce~ Xt (INI^ + a 2 (0) + 6 2 (0)) . (4.27) 



5. Proof of Theorem HB 

5.1. Preliminary definitions. We first establish some definitions and notation. 
For 7 (x) € M 2 with components a(x) and (3{x) denote | 7 (ai)| = |a:(x)| + |/?(x)|, and 

iiTiioo = sup | 7 (x)i, n 7 || zl = [mm. (5.1) 

xe[o,i] Jo 

In what follows, for a time-varying vector 7(x, t), we denote | 7 | — | 7 (x, i)| and ||7|| = 
|| 7 (-,i)|| to simplify our notation. For a 2 x 2 matrix M, denote: 

\M\ = max{|M 7 |; 7 G M 2 , | 7 | = 1}. (5.2) 

For the kernel matrices K{x, £) and L(x,£) denote 



sup \K{x,£)\. 

(x,e)£r 



(5.3) 



For 7 € -ff 2 ([0, 1]), recall the following well-known inequalities, that will be used later: 



||7l|oo<C 3 [|| 7 || i2 + || 7;c || L 2]<C4|| 7 || ffl , 

||7x||oo < C 5 [||7 x ||i2 + HtkxIU 2 ] < CellTllja-2. 



(5.4) 
(5.5) 
(5.6) 



Define the following linear functionals, the first two of which are, respectively, the 
inverse and direct transformations (I3.23[) and (|3. 381) : 



1C[j](x 
£[ 7 ](x 

}Ci[i\{x 

/C 2 [ 7 ](x 
A [71(1 

£u[i\(x 



Jo 

j{x,t) + { L(x,0^,t)d^, 
Jo 

-K(x, x) 7 (x, t)+ [ K^x, £h(£, t)d£, 
Jo 



= -K{x, x) 7 (x, t) - / K x (x, Ol(t t)d£, 
Jo 

= L(x,x) 7 (x,f) + / L x (x,0l(Z,t)d£, 
Jo 

= (L x (x,x) + L^(x,x))^(x,t) + L x (x,x)j(x,t) 



+ / L xx (x,Ol(^t)dC 

'0 



(5.7) 
(5.8) 
(5.9) 
(5.10) 
(5.11) 

(5.12) 
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For simplicity, in what follows we drop writing the x dependence in functionals and 
the t dependence in the variables. 

Using (pT7)) and |[5T5]) . we define Pi [7] and ^[7] as: 

F 1 =K NL {L[ 1 ],x), F 2 = f NL (£fr],x). (5.13) 

To prove Theorem 14.11 we notice that if we apply the (invertible) backstepping 
transformation (|3.23p to the nonlinear system (|4.8p we obtain the following trans- 
formed system: 

= 7t - S(a;)7x + Ajvl(w, 2;)^ + Inl(w, x) 

+ [ K(x,0(A NL (w,Owx(0-fN L (w,0)dt;, (5.14) 



and using the inverse transformation (13.38)) the equation can be expressed fully in 
terms of 7 as: 

7t - E(x) 7x + F 3 [-y,<y x ] + F 4 [ 7 ] = 0, (5.15) 

where the functionals P3 and P4 are defined as 

P a =/C[Pi[7hJ, (5-16) 
F4=/C[Fi[ 7 ]£i[7]+F 2 [ 7 ]]. (5.17) 

The boundary conditions are 

o(0, t) = 9/3(0, t) + G NL (/3(0,t)) , /9(l,i) = a(t) + &(t). (5.18) 

By the assumptions on the coefficients and applying Theorem IA.21 the direct and 
inverse transformations f|3 . 23[) and (|3.38[) have kernels that are C 2 (T) functions. Dif- 
ferentiating twice with respect to x in these transformations, it can be shown that the 
H 2 norm of 7 is equivalent to the H 2 norm of z (see for instance [33J ) . Thus, if we 
show H 2 local stability of the origin for (|5.15p - (|5.18|) . the same holds for z. 

We proceed by analyzing (using a Lyapunov function) the growth of || 7 11^,2, || 7*11^2 
and ||7tt||L 2 - Relating these norms with ||7||^2, we then prove H 2 local stability for 
7- 

5.2. Analyzing the growth of ||7||l 2 - Define 

V 1 = [ -i T {x,t)D{x)-i{x,t)dx, (5.19) 



for D(x) as in (|3.9p . Proceeding analogously to (|3.1ip - (|3.14p , we get some extra 
nonlinear terms: 

Vi = - l T (x, t) (D(x)Z(x)) x j(x, t)dx + [ 7 T (x, t)D{x)Y,{x)l{x, t)] I 



7 T (x, t)D(x) (P 3 [7,7x] + PiW) dx. (5.20) 







Let us analyze first the last term: 
-1 

' <Ki / |7|(|P 3 [7,7JI + I^[7]|)^. 



7 T (x,i)i?(x) (P3[7, 7 J + P 4 [7])^ 







(5.21) 
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Applying Lemma IB.2I (see the Appendix) , we obtain that there exists a 8\ , such 
that for ij-ylloo < 5 X , 

H|i ? 3[7,7x]|da : <^a||7«lloo||7lli», (5-22) 
l7l|f4[7]|^<^3ll7l|oo||7l| 2 L2, (5.23) 

1 /2 

and using inequality (|5.5p and noting that ||7||l 2 < K^V X , we obtain 

{FafaixWdx < KsWkWvoVi, (5.24) 

| 7 ||F 4 [ 7 pz < ^6||7x||oo^i +K 7 V 1 3/2 . (5.25) 



Now, 

[7 T (x, t)D(x)S(a:)7(a;, t)] ] = Ba 2 (t)e^ - Aa 2 (l, t)e _ " - 5/3 2 (0, t) 

+A( 9 /3(0,i) + G AfL (/?(0,t))) 2 , (5.26) 
and for HtH^ < 5i, \G NL (0(O,t))\ < K 8 \/3(0, t)\, and A > 0, we obtain 

[j T (x, t)£>(a;)E(a:)7(a 1 1)] J < -Aa 2 (l, t) e -" + + X 8 ) 2 - 5)/3 2 (0, t) 

+Ba 2 (t)e fi (5.27) 

Thus, choosing B = (\q\+ Kg,) 2 A + A 2 and A and /i as in the proof of Proposition 13. 11 
we obtain the following proposition: 

PROPOSITION 5.1. There exists S± such that if ||7||oo < #1 then 

Vx < -A1V1 - A 2 (a 2 (l, t) + /3 2 (0,i)) + C X V^ 2 + C 2 hx\\ooV 1 

+C 3 {a 2 (t) + b 2 (t)), (5.28) 

where X\, A2, Ci, C2 a?jd C3 are positive constants. 

5.3. Analyzing the growth of 1 1 Tt II - Define 77 = j t - Notice that the norms 
of r\ and j x are related (see Lemma lB.6l in the Appendix). Taking a partial derivative 
in t in f|5 . 14[) we obtain an equation for rj as follows: 

r) t + (Fi[ 7 ] - E(ar)) ^ + F 6 [7,7x,»?] + F 6 [ 7 ,»7] = 0, (5.29) 
where F$ and Fq are defined as 

F 5 = K x \F X YlH + f X K(x, £)F 12 [ 7) jMOdC + K(x, 0)A NL ( 7 (0), 0) r?(0) 
Jo 

+/C[F 11 [ 7 ,r 7 ] 7c ], (5.30) 
F 6 = K [FnfaTilC [7]] + AC [F x [ 7 ]£, [ry]] + /C[F 21 [ 7 , 77]], (5.31) 



where 



a* 

F 11 = ^^(C[ 7 ],x)C[r ] }, (5.32) 

F 12 = ^ (C[ 7 ],x) ( 7x + C X H) + (£[ 7 ],z) , (5.33) 

F 21 = ^(£[ 7 ] )a; )LM. (5.34) 
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The boundary conditions for r\ = \r\\ rj 2 ] T are 

T?i(0,i) = qm(0,t) + G' NL (f3{0,t)) m {0,t), r) 2 (l,t) = -d ia (t) - d 2 b(t). (5.35) 
To find a Lyapunov function for rj, we use the next lemma: 

Lemma 5.2. There exists S > such that, for H7H00 < S, there exists a symmetric 
matrix i?[7] > verifying the identity: 

R[ 7 ] (£(*) - Fx [7]) - (E(x) - F [ 7 ]) T J2[ 7 ] = 0, (5.36) 
and the following bounds: 

R[i\{x) < ci + oalMU, (5-37) 

I ((i?[ 7 ] - £>(*)) £(*)), | < c 2 || 7 ||oo (1 + || 7 *|U , (5.38) 

\(Rh}) t \<cs(\v\ + H\^), (5.39) 

where Ci,C2,C3 are positive constants. 

Proof. We explicitly construct i?[ 7 ] as 

R[j] = D{x) + 6[ 7 ], (5.40) 

with 



e[ 7 ] 

where ^[7] is defined as: 





^[7] 



(5.41) 



/ r i _ gu(g) (Fi[ 7 ]) 12 -^2 2 (^) (gtfrDgi 42 ^ 
WJ £2 (a:)+ ei (x) + (F 1 [7]) 11 -(F 1 [7]) 22 ' 1 ■ J 

where (Ft [7])^ denotes the coefficient in row i and column j in the matrix Fl[ 7 ]. 
Identity f|5 . 36[) follows by using the construction of i?[ 7 ](x) in (|5.40p - (|5.42p . and the 
fact that D{x) and £(x) are diagonal and commute. To ensure that the denominator 
of (|5.42| is different from zero, denote K\ = minimi] (ei(x) + £2(2^)) > 0. Applying 
(|B.13|) from Lemma fB.2| there exists Si for which, if || 7 ||oo < Si, one gets: 



e 2 (x) + ei(x) + (Fi[ 7 ]) 22 - (F 1 [ 7 ]) u > K x - K 2 \\ 7 \U (5.43) 

2K 2 ■ 



thus if I j 7 1 j oo < minj^!,^} with 5 2 = -M^-, we obtain 



e 2 ( 2 ;)+e 1 (x) + (F 1 [7]) 22 -(F 1 [7]) 11 > ^1, (5.44) 

thus ip[y] is well-defined. Applying again (|B.13[) in the numerator of (15.42)) to bound 
(|5.4ip . we obtain: 

||6[ 7 ]|| <if 3 || 7 ||oo, (5.45) 

and noting K4 = ||-D||oo, we obtain directly the bound (|5.37l) . and by choosing H7II00 < 
min{<5i, S 2 , S3}, with 5 3 = we show R[y] > 0. 
Inequality (|5.38p is equivalent to showing: 

|(e[7]S(x))J< C2 || 7 ||oo(l + ||7,||oo). (5.46) 
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We first use (|5.45[) to bound |0[7](a;)I] a .(a;)|, and for |9 x [7](:e)E(;e)| we take a deriva- 
tive in x in (|5.42p . use the bound (|5.44[) and use the fact that -§^.Fi['j\ = -Fi2[7, 7J 
and using Lemma |B.3[ there exists 84 such that if H7II00 < 84, 

\Fi2bn x ]\<Kx(H\oo + h x U)- (5-47) 

To show (|5.39p we use J^Fi[7] = | i^i 1 [t, 77] and apply Lemma 133.31 Setting 5 = 
min{(5i, S 2 , S3, 64} the lemma follows. □ 
Define: 

V 2 = f r 1 T (x,t)R[ 1 ](x)r ] {x,t)dx. (5.48) 
Jo 

Computing V2, applying Lemma IS~!?I and integrating by parts, we find 
V 2 = - f\ T (x,t) (R[j\ (E(z) - F 1 [j])) x i](x,t)dx 



[ V T (x, t)R[y] (x) (E(z) - F 1 [7] (a;)) V (x, t)] X J Q + I if{x, t) (R{ 7 }\ V (x, t)dx 



-2/ i 1 1 {x,t)R[ 1 ]F b [ 1 , lx ,r ] ,r lx ,)dx^2 r) 1 (x,t)R[i\F 9 tr,rj\dx. (5.49) 



The first three terms of (|5.49[) are analyzed using Lemma \b. 2 1 Thus, there exists 8\ 
such that, for I^Hoo < 5, we find, for the first term: 

1 

V T (x, t) (i?[ 7 ] (£(x) - Fx [ 7 ])) x V (x, t)dx 

< -\ x v 2 + KtWnWl, (|| 7 || 00 H~ 1 1 'fx 1 1 00 

) . (5.50) 

The second term of (|5.49l) is bounded using the boundary conditions, (|3.1ip - (|3.14p . 
and Lemma 15.21 as: 



[ V T (x,t)R{~f}(x) (ECaO-iMlK^M^Clo 



< -A 2 (r)l(l,t) + 4(0,1)) +K2W7W0C ( V 2 2 (0,t)+ V 2 1 (l,t)) 

+A' 3 (l + ||7lloo)(a(t) 2 + 6(i) 2 )- (5-51) 
Finally, we bound the third term of (|5.49|) applying Lemma [S~2l as follows: 



V T (x,t)(R[y]) t r 1 (x,t)dx<K 3 / \ v \ 2 (\ v \ + \\r,\\ L i) dx < K4\\v\\U\vU- (5.52) 

Jo 

Applying Lemmas 15.21 and IB. 31 to the last terms of (|5.49[) , we get, for ||7||oo < 



<k 5 / inWF^njldx 



r) T (x, t)R[j]F 5 [7, j x , 77, r] x }dx 
< KeWvWh (hlloo + H7.II00) + K 7 \\v\\M0,t)M0,t)\, (5.53) 

and 



1.1 



7 ? T (a;,i)i?[7]F 6 [7,7 y ]da; 



1 



<K S \v\\Feh,v]\dx < KsWvWUhlU (5.54) 
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Thus, it is clear that by choosing ||'y || oo small enough, using Lemma lB.6l to bound 
IITsIIoo by 1 1 77 1 1 , and noting ||??||l2 < K W V 1 / 2 ^ we obtain the following proposition: 
PROPOSITION 5.3. There exists 6 2 such that if \\j\\oo < S 2 

V 2 < -\ 3 V 2 -\ i ( V 2 1 (l,t)+r ] 2 2 (0,t))+K 1 V 2 \\i 1 \\ oo + K 2 b(t) 2 , (5.55) 

for X 2 , A3 , K\ , K 2 positive constants. 

5.4. Analyzing the growth of ||7tt||£2. Define 9 = r\ t . Notice that the norms 



of 9 and r] x are related (see Lemma lB.8l in the Appendix). Taking a partial derivative 
in t in (|5.29| we obtain an equation for 9: 

t + (Ft [7] - E(a;)) 9 X + F 7 [ 7 , 7*, »?, Vx,&] + F s [7, V, 0] = 0, (5.56) 
where F-j and Fs are defined as 

F 7 = JC 1 [F u [ 1 ,r ] ]r,}+ I K(xX)F l2 [ lllx ]9(^)di + lC l [F l [ 1 n 
Jo 

+ f K(x, 0F 14 [ 7 , 7*, V, vMOdt + K(x, 0)^^ ( 7 (0), 0) r?(0)r/(0) 
Jo 07 

+K(x, 0)Anl (7(0), 0) 9(0) + JC [Fu [7, vM + K. [F 13 [j, V, %*] , (5-57) 

F 8 = 2JC [F n [7, rj\C x + JC [F x [y\C x [9]] + K [F 13 [ 7 , 77, 9]C X [7]] 

+/C[F 22 [ 7 ,77,0]], (5.58) 

where 

F i3 - ^ (C[~<),x) L[ri\C[rj\ + ^ (£[7],x)£[0], (5.59) 

r) 2 A r)A 

^ = (£[7],*) (7, + c x [i\) + -f± (Cfrlx) ( Vx + cm) 



d 2 A 



T-Q^L(<Ch],x)£[r,}, (5.60) 

F22 = (Cb\,x)C[rj\C[rj\ + ^ (£[7], 1) £[0]. (5.61) 

The boundary conditions for # = [9\ 9 2 ] T are 

0i(O, t) = q9 2 (0, t) + G' NL (/3(0, t))9 2 (0, t) + G" NL (/3(0, t))r,l (0, t), (5.62) 

9 2 (l 1 t) = d\a(t) +djb(t). (5.63) 

Since (|5.56[) has the same structure as (|5.29[) . we define: 

V 3 = [ 9 T (x,t)R[ 1 ](x)9(x,t)dx, (5.64) 



where R[y](x) was defined in Lemma [ 

Computing V3, and proceeding exactly as in (|5.49l) . we find: 



1 

Tf 



V 3 = - 9 1 (x,t)(R[ 1 ](Z(x)-F 1 { 1 })) x 9(x,t)dx 
Jo 

+ [9 T (x, t)R[ 7 ](x) (£(*) - Fx Ij](x)) 9(x, t)] ^ + [ T (x, t) (R[ 7 ]) t 6(x, t)dx 

Jo 

9 T (x,t)R[ 1 }F 7 [ 1 , lx , V ,r lx ,9}dx-2 9 T (x,t)R[- f }F 6 [ 1 ,n,9]dx. (5.65) 
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The first three terms of (|5.65[) are analyzed as in (|5.50l) - (|5.52l) : 

V 3 < -X1V3 + KiWOfv (hlloo + ||T«I1oo) + MtIIoc - A 2 ) + 9 2 (0,t)) 



+2 



9 T (x,t)R[ 1 ]F 6 [ 1 ,r,,e]dx 



(x,t)R['y]F 7 ['y,j x ,r),r) x ,9]dx 



-K 3 \\nl4v\\oo + Mrfi(0,t) + (1 + ||7l|oo)(a 2 (t) + b 2 {t))). 



(5.66) 



Finally, applying Lemmas 15.21 and IB. 41 in the last two terms of (|5.66p , there exists a 
S, such that for IMIoo < 5, 



(x,t)R[y]F 7 [<y,j x ,ri,r] x ,6}dx 



leWFri^j^^d^dx 



< K e \\e\\ 2 L2 (hlloo + hxlloo) + K 7 \\6\\ L 4 V \\l 2 + tfe||0|MWMM|ao + K 9 \\9\\ L2 
+K 1Q \\9\\ L > {\M\l4v\L + \v(0,t)\ 2 + | 7 (O,t)||0(O,t)|) , (5.67) 

and 



9 T (x,t)R[ 7 }F 8 [ 1 ,r ] ,e}dx 



< K 
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\e\ \f 8 h,v,o]\dx 



< ^■nll^lli^llTlloo + ^lall^lU-H^IU^II^Iloo + ^all^lU-ll^lli^ + ^i 3 ||»?||i 2 ||^IU-.(5.68) 

Thus, by choosing H7II00 and ||??||oo small enough to apply Lemma |B. 81 we finally 
obtain the following proposition: 

Proposition 5.4. There exists S3 such that */ UtIIoo + ll^llcn < ^3 then 

v 3 < -hv 3 - a 6 (ef(i,t) + e 2 2 (o,t)) + k^v^ 2 + k 2 v 2 v 3 1/2 

+K 3 V 3 3/2 + if 4 h||oo% 2 (0,t) + K b {a 2 + b 2 ), (5.69) 

where A5, Xq, K\, K 2 , K 3l K4, K§ are positive constants. 

5.5. Proof of H 2 stability of 7. Defining W = Vi + V2 + V3, and combining 
Propositions 15.11 15.31 and !5.4[ there exists 6 such that if ||7||oo + IMloo < 8 



W < -XiW + GW 3/2 + C 2 (a 2 + b 2 ), 



(5.70) 



for X,C\,C 2 > 0. To compensate the last term, we augment this Lyapunov function 
and define 5 = W + f (^ + £). Then, 



5 < -XiW + CW 3/2 + (C 2 - c)(a 2 + & 2 ), 
and choosing c > C 2 , one obtains 

5 < -X 2 S + dS 3/2 , 



(5.71) 



(5.72) 



for some positive A 2 . Following [2] and noting H7H00 + || r /||oo < C 2 S. then for suffi- 
ciently small 5(0), it follows that S(t) — > exponentially. 

Given that W (by Proposition IB.5|) is equivalent to the H 2 norm of 7 when 
II 7 II 00 + IMloo is sufficiently small, and since by construction 70 verifies the required 
second-order compatibility conditions, there exists 6 > and c > such that if 
||7o||iP < s , tnen: 



||, + a {tf + b(t) 2 < ce" At (||7o|j 2 H2 + a(0) 2 + b(0) 2 ) 



(5.73) 
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Since, as we argued, for small enough ||^||//2 the H 2 norms of z and 7 are equivalent, 
this proves Theorem 14. II 

Remark 2. The proof has been carried out for the case q 7^ 0. If q = 0, we have 
to modify the target system following Section 13.51 and this implies the appearance of 
a linear boundary term (a coefficient times /3(0, t)) in the 7 system; similarly, in the 77 
and 9 systems, 772 (0, and #2(0,^) terms will appear. These terms can be controlled 
using the same Lyapunov function by following the strategy outlined in Section [ 



6. Concluding remarks. We have solved the problem of full-state boundary 
stabilization for a 2 x 2 system of first-order hyperbolic quasilinear PDEs with actu- 
ation on only one boundary. We have shown, using a strict Lyapunov function, H 2 
local exponential stability of the state. It is possible to extend this result to design an 
observer, as shown in (36], and combining both results one obtains an output- feedback 
controller with similar properties (see |37|). 

It would be of interest to extend the method to n x n systems. For instance, a 
3x3 first-order hyperbolic system of interest is the Saint- Venant-Exner system, which 
models open channels with a moving sediment bed (7| ; the extension is shown (for the 
linear case) in [9]. While extending the Lyapunov analysis ton xn systems has been 
done [3j, considerable extra effort is required to extend backstepping to a general nxn 
system, even in the linear case. In general, the method needs n 2 kernels resulting in 
a n 2 x n 2 system of coupled first-order hyperbolic equations, whose well-posedness 
depends critically on the exact choice of the transformation and target system. The 
extension has been shown possible, for the linear case, if the system has n positive 
and one negative transport speeds, with actuation only on the state corresponding to 
the negative velocity |10] , 

Appendix A. Well-posedness of the kernel equations. We show well- 
posedness of the following hyperbolic 4x4 system, which is generic enough to contain 
all the kernel equation systems that appear in the paper: 



e 1 {x)Fl + ei {i)Fl = 91 (x, £) + J2 C u (x, £)F* (x,0, (A.l) 

4 

e^F 2 - e 2 (0F 2 = g 2 (x,0 + J^CwfaO^foO. (A.2) 

4 

e 2 (x)F 3 - ei (0if - 93(x,0 +^C 3i (x,0F i (x,$, (A.3) 

i=i 
4 

e 2 (x)F* + e 2 (0F* = 9i (x,0 + £ C^foO^.O, (A.4) 



evolving in the domain T = {(x, £) : < £ < x < 1}, with boundary conditions: 

F 1 (x, 0) = hi(x) + qi(x)F 2 (x, 0) + q 2 (x)F 3 (x, 0), (A.5) 
F 2 (x,x) = h 2 (x), F 3 {x,x) = h 3 (x), (A.6) 
F A (x, 0) = h 4 (x) + q 3 (x)F 2 (x, 0) + q 4 (x)F 3 {x, 0). (A.7) 



This type of system has been called "generalized Goursat problem" by some au- 
thors |19j. However the boundaries of the domain T are characteristic for (|A. 1|) 
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and (|A.4[) . thus the general results derived in [TS| cannot be applied. The following 
theorems discusses existence, uniqueness and smoothness of solutions to the equations. 

Theorem A.l. Consider the hyperbolic system HA. l\) - [A~7fy . Under the as- 
sumptions qi,h{ G C([0, 1]), gi,Cji G C(T), i,j — 1,2,3,4 and e\,ti G C([0, 1]) with 
e i( x )i e 2(x) > 0, there exists a unique C(T) solution F l , i = 1,2,3,4. 

Theorem A. 2. Consider the hyperbolic system liA. l]) - fA. 7j) . Under the assump- 
tions of Theorem ] A. 1[ and the additional assumptions €i,qi,hi G C N ([0, 1]), gi,Cji G 
C N (T), there exists a unique C N (T) solution F\ i = 1, 2, 3, 4. 

Next we prove the theorems; the proof is based on transforming the equations into 
integral equations and then solving them using a successive approximation method. 

A.l. Transformation to integral equations. The equations can be trans- 
formed into integral equations by the method of characteristics. For that, it is neces- 
sary to define: 

f x 1 f x 1 

M x ) = / —^-dz, M x ) = / —j-^-dz, (A.8) 
Jo e n z ) Jo e 2V z ) 

and 03 (x) = 4>i(x) + 02(x). Note that all the functions are monotonically increasing 
and thus invertible, due to positivity of the e coefficients. Under the assumptions of 
Theorem |A~T| it also holds that 0*, 0" 1 G C 1 ([0, 1]). 

Define, for (x,£) G T, the characteristic lines along which (jA.lj) - ()A.4j) evolve: 



Xl (x,Z,s) = fc 1 (Mx) ~ MO + s) , (A.9) 

^(x^,s)=cf>^(s), (A.10) 

x 2 (x,£,s) = M (0! (03 1 (Mx) + M0))+*) , (ATI) 

s) = 0J 1 (0 2 (03 1 (0!(X) + MO)) ~ *) , (A.12) 

X 3 (X,£, s) = 2 - X (02 (03' 1 (02 (X) + MO)) + S) , (A.13) 

& s) = M (0! (03' 1 (0 2 (x) + MO)) ~ s) , (A.14) 

x A (x, f , s) = 0^ (Mx) - MO + s ) > ( A -15) 

U(x^,s)=^ 1 (s), (A.16) 

where the argument s that parameterizes xi and £j belongs to the interval [0, sf], 
with sf defined as 

s((x,0 = M0, (A.17) 

4{x,£) =0i(x)-0 1 (0 3 ' 1 (0 1 ( a; )+02(O)), (A.18) 

sf(a;,0 = 02(x)- 2 (0 3 - 1 (02(x)+0i(O)), (A.19) 

«f(s,0 = fc(0- (A-20) 
The following holds 



Lemma A. 3. If(x,£) G T ands G [0, sf], itholds that (xi(x, £, s), £i(x, £, s)) G T, 
for i — 1,...,4. Also, under the assumptions of Theorem \A.1[ X4, and sf are 
continuous in their domains of definition since they are defined as compositions of 
continuous functions. Moreover, the following inequalities are verified 



Xi(x,£,s) <x, i = 1,...,4, 



(A.21) 

(A.22) 
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Using these definitions, (|A.l[) - (jA.4[) are integrated to: 

Fi(x,0 = F j (x j (x,Z,0)^ j (x,Z,0)) + G j (x,Z)+I j [F](x,O, (A.23) 
where we have denoted 



F = 



F 1 
F' 2 
F 3 
F i 



G j (x,£) = 







9j (xj(x,Z,s),Zj(x,Z,s))ds 1 (A.24) 



Cji (xj (x, Z, s), Zj (x, Z, s)) F l (xj (x, Z, s), (x, f , s)) ds,(A.25) 





for j = 1,2,3,4. Substituting the boundary conditions (|A.5|) (|A.7[) and expressing 
the terms in (|A.5|) and (|A.7j) containing F 2 (x, 0) and F 3 (x, 0) in terms of the solution 
(|A.25[) we get four integral equations which have the following structure: 

F* (x, = Hj (x,0 + G (x, + <p d (x, + Qj [F] (x, Z) + I [F] (x, Z), (A.26) 

where Hj(x,Z) = hj(xj(x,Z,0)), (fj(x,Z) is has the values ^ = f3 — and 

<Pi = qi{xt(x, Z, O^^x, £ 0), 0) + 92(11 (x, Z, 0))H 3 ( Xl (x, Z, 0), 0) 

+q 1 (x 1 (x, Z, 0))G 2 (xi(x, C, 0), 0) + teMz, 6 0))G 3 (xi(x, £, 0), 0), (A.27) 

+q 3 (xi(x, Z, 0))G 2 (x 4 (x, Z, 0), 0) + qi {x A {x, Z, 0))G 3 (x 4 (x, f , 0), 0). (A.28) 
and the values of the <3j[-F](x, £) are Q 2 — Q3 = and 

Qi [F] = qi ( Xl (x, Z,0))I 2 [F] ( Xl (x, Z, 0), 0) + g 2 (si (x, £, 0))/ 3 [F] (x x (x, 0), 0), (A.29) 
Q 4 [F] = 93 (^4 (x, Z, 0))/ 2 [F](x 4 (x, Z, 0), 0) + g 4 (x 4 (x, Z, 0))/ 3 [F](x 4 (x, Z, 0), 0). (A.30) 

In this form, the equations are amenable to be solved using the successive ap- 
proximation method. This is explained next. 

A. 2. Solution of the integral equation via a successive approximation 
series. The successive approximation method can be used to solve the integral equa- 
tions. Define first the following functional acting on F: 



$ j [F}(x,Z) = Q j [F}(x,Z) + I j [F](x,0, 



and the vectors: 



H x + Gi + pi 

H 2 + G 2 + lf2 

H 3 + G 3 + <p 3 

H4 + G4 + (f4 



$[F] 



$i[F] 
$ 2 [F] 
$ 3 [F] 



(A.31) 



(A.32) 



Define then 



F°(x, Z) = <p{x, Z),F n (x, = flF"- 1 } (x, 0- 



(A.33) 
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Finally define for n > 1 the increment AF™ = F n — F n 1 , with AF° = <p by defini- 
tion. It is easy to see that, since $ is a linear functional, the equation AF n (x,0 = 
<f>[AF n - l ](x,ti) holds. 

If linin^oo F n (x, £) exists, then F = linin^^ F n (x, £) is a solution of the integral 
equations (and thus solves the original hyperbolic system). Using the definition of 
AF n , it follows that if YlnLo AF n (x, £) converges, then 

oo 

F(x,0 = Y, AF "( x >0- (A.34) 

n=0 

A. 3. Proof of convergence of the successive approximation series. First, 
define: 

4> = , m ax {\<fii(x, 01} , Cji = max \Cji(x, £)| , K e = max j 

(*>£)er (is.Oer (*,f)eT [ ei(a;) e 2 (a;) 

i=l,2,3,4 

g< = max \ qi (x)\ , C = ( 1 + 9i J ( E E <^ I ' ( A ' 35 ) 

Next, we prove the following two lemmas: 

Lemma A. 4. For i = 1,2, 3,4, n > 1, (x, £) e T, and sf (x, £), Xi(x, £, s) defined 
as in iA.9f) - fA.20\) . it follows that 



'(*,«) 



(i 



x?(x, £ s)ds < K t —^. (A.36) 



Proof. We show the result for i = 1, 2. It follows for i = 3, 4 by switching e! and 
(>i, respectively, for e 2 and <j) 2 . For i = 1 we can write: 

<M0 rMO 

x?(x,Z,s)ds= [ft 1 (fa(x) - MO + *))] da. (A.37) 

Jo 

To prove the inequality, change the variable of integration to z — 0f 1 (4>i(x) — <pi(0 + s) 
Then, taking into account 

£ = Ts ^ {Mx) " M0 + s)] = ^) = £l(z) ' (A ' 38) 

the integral can be bounded as follows: 

r<M0 



/■x (-a; T n+1 

/ z n /e 1 {z)dz < K e z n dz = K € -. (A.39) 



i<j>l'(Mx)-Mi)) 
For i = 2 the integral can be written as: 



o 



0i(a : )-0 1 (0- 1 (0 1 (x)+0 2 (C))) 

[M 1 (fa (C 1 (0i(») + &(£))) + s )f (A.40) 
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As before, change the variable of integration to z = 4>\ 1 (01 (03 1 (4>i{ x ) + 02(C))) + si- 
Then one has that ^| = -^rjpj — ^i(z), thus the integral can be bounded as follows: 



<p 1 (x)-<t> 1 (4>3 1 (<t>i(x)+Mt))) 

[C 1 (0i (03 1 (M*) + MO))+s)] n ds 



0- 1 (0 1 (x)+0 2 (^)) 



/•z T n+1 

z n /e 1 (z)dz < K e / z"dz = K e 

Jo n + l 



which concludes the proof. □ 

Lemma A. 5. For i = 1, 2, 3, 4, n > 1 and (x, £) 6 T, assume that 



(A.41) 



|Ai7(z,OI<0- 



(A.42) 



t/ien it follows that \<$>i{AF n ){x, £)| < 0^— (^tTyr— ■ 

Proof. We show it for i = 1,2; the structure of the equations is the same for 
i = 3,4. For i = 2: 



|$ 2 [AF"](x,OI = \I 2 [AF n ](x,0\ <J2°* 



i=l 



' (*,€) 



\AF? (x 2 (x,S,3),b(x,£,8))ds\ 



n 



2 1 



1=1 



' (*>«) 



g-n+l Qn+l g.n+1 

(n + 1)! 



, (A.43) 



where Lemma lA.41 has been applied. Similarly, for i = 1: 



^[AF n ](x,0 < IQxtAJ^K^OI + |/i[AF"]0r,OI 

4 r s|'(xi(s,e,0),0) 



< 910- 



a;2(a;i(a;,^,0),0,s)ds 



sf (zi(x,^,0),0) 



a#(aii(a:,f,0),0,s)ds 



+0- 



4 



x"(x, £, s)ds 



< 910 ; V 2 j — h 920 : > C 3i 

n. n+1 n! ' 



n! ^ n+l ~ 



71 +1 71! 



(n + 1) 



n+l 
(A.44) 



since Xi(x,£, 0) < x. Thus the lemma is proved. □ 
Next we show that (|A.34|) converges. 

Proposition A. 6. For AF?(x,£), i = 1,2,3,4, one has that 



< (be 



CK e x 



(A.45) 



n=0 
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Proof. The result follows if we show that \AF?{x,£)\ < <t> „f - ■ We prove the 
bound by induction. For n = 0, the result follows from (IA.33|) . Assume that the 
bound is correct for all i in AF n (x,£), Then, we get for AJ 7 ™ +1 (x,^) that 



_ (~m+l Z(n+1 T n+1 

\AF?+\x,t)\ = |$,[AF"](x,OI < 4> , (A.46) 

1 1 (n + 1)! 



where we have used Lemma lA. 51 Thus the proposition follows. □ 



From Proposition IA.6I we conclude that the successive approximation series is 
bounded and converges uniformly. Thus, a bounded solution to Equations (|A. 1|) 
(|A.7p exists. This proves the existence part of Theorem lA.il 



To prove uniqueness, let us denote by F(x,^) and F'(x,^) two different solutions 
to (|A~T1l- ([A~7|) . Defining F(x,£) = F(x,£) - F'(x,£). By linearity of (|AlTl - fA~7|) . 
F(x,£) also verifies (|A.1 |) — (|A.7 |) ■ with h t = for all i. Then cj> = for Ffag), and 
Proposition I A. 61 we conclude F(x,£) = 0, which implies that F(x,£) = F'(x,£). 



To prove that the solution is continuous, note that since (|A.34[) converges uni- 
formly, one only needs to prove continuity of each term. First, AFq = ipi G C{T) 
since the ifi are defined as a sum of compositions of continuous functions. Similarly, 
since AF n is defined as the integral (with continuous limits) of continuous functions 
times the previous AF n ^i composed with continuous functions, by induction it can 



be shown that AF n € C(T). Thus F G C(T) and Theorem lATl is proved. 



A. 4. Smoothness of solutions. Next we sketch the proof of Theorem IA.2I 
We only consider N = 1; for N > 1 the result can be proven by induction. Denote 

dFj 
3€ 



Gi = Qp-(x,t;) and Hi = By differentiating with respect to x and £ in 
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(|A.ip - (IA.4|) . we find two uncoupled 4x4 hyperbolic systems (for G, ; and for Hi) 



e l (x)Gl + e 1 (OG\ = -e'^F 1 + ^(x,0 + £ ^-{x^M 

i—1 



(A.47) 



i=l 



i—1 



(A.48) 



i=l 



i—1 



e 2 (x)G^ + e 2 (0^ 



^C 3j (x,e)G l (x,0, 



(A.49) 



i=i 



i=l 



■J2C u (x^)G i (x,0, 



(A.50) 



i=l ' 



(A.51) 



i=l 



i=l 



(A.52) 



e 2 (x)H^-e 1 (0Hl 



-^)F 3 + ^(x,0+T,^T^0F\x,0 

i=l ' 



e 2 (x)Ht + e 2 (0Hl 



(A.53) 



i=l 



i=l ' 



(A.54) 
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Now, differentiating the boundary conditions in (|A.6[) it is found that: 

G 2 (x,x) + H 2 (x,x) = h' 2 (x), G 3 (x,x) + ff 3 (x,x) = h' 3 (x), (A.55) 
and setting x = £ in (|A.2|) - (|A.3j ) 

4 

ei(x)G 2 (x,x) - e 2 (x)ff 2 (x, x) = g 2 (x,x) + ^ G 2l (x, x)F l (x, x), (A.56) 

i=l 
4 

e 2 (x)G 3 (x, x) - ei(£)# 3 (>, i) = <7 3 (x, x) + £ G 3l (x, x^x, x), (A.57) 

i=l 

we can find a set of boundary conditions for G J and fP, j = 2,3, at the boundary 
x = i: 

G 2( x x = e 2 (x)h' 2 (x) + g 2 (x,x) +J2t=i C 2i(x,x)F i (x,x) ^ 

e 2 (x)+ei(x) 

^3/ \ _ ei(x)h' 3 (x)+g 3 (x,x) +Y l t=i C 3i(x,x)F i (x,x) 

& (X ' X) ~ e 2 (x)+ £l (x) ' (A ' 59j 

„2, s _ ei(x)h' 2 (x)- g 2 (x,x) - J2t=i C 2i{x,x)F i (x,x) 

a [X,X) — — — ■ — , (^A.DUJ 

e 2 (x) + ei{x) 

tt3( n _ e 2 (x)h' 3 (x)- g 3 (x,x) - G 3i (x, x)F'(x, x) 

e 2 (x) + ei(x) 

Similarly, differentiating boundary conditions (|A.5[) and (|A.7[) we find boundary 
conditions for G 1 and G 4 at £ = 0: 

G l (x,0) = h[(x) + q' 1 (x)F 2 (x,0) + q' 2 (x)F 3 (x,0) + qi (x)G 2 {x,0) 

+g 2 (x)G 3 (x,0), (A.62) 

G 4 (x, 0) = fti(a?) + g 3 (x)F 2 (x, 0) + ^(x)F 3 (x, 0) + g 3 (x)G 2 (x, 0) 

+q 4 (x)G 3 (x,Q), (A.63) 

and setting £ = in ([A~T|) - (|X4|) we can also find two sets of boundary conditions for 
iP, j = 1,4, at the boundary £ = 0: 

!(x, o) = giMRS^foo^M) _ fig) ( K{x) + q[{x)F 2 (x , 0) 



£l (0) ei (0) 

giQ) 

ei(x) 



+g 2 (x)F 3 (x, 0) + U(0)# 2 (x, 0) + 52 (x, 0) + ]T G 2j (x, 0)F*(x, 0) 



12 (x) 

e 2 (x) 



ei (0) J ff 3 (x,0)+ 5 3^,0) + ^G 3j (x,0)F l (x,0) , (A.64) 



H\x, 0) = g 4 ^0)+Et 1 ^(^0)^(x,0) _ e^) / + 

£2(0) e 2 (0) V 

+^(^)^ 3 (^, 0) + (e 2 (0)H 2 (x, 0) + 52 (x, 0) + J2 C *( x > °) Fl ( x > °) 



?4(x) 

e 2 (x) 



ei(x) 

e 1 (0) J ff 3 (x,0)+ 53 (a ; ,0) + ^G3 l (x,0)F l (x,0) j j . (A.65) 
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Thus, both the G l 's and H l, s verify equations formally equivalent to (|A.ip - (|A.7[) , 
with derivatives of the old equations' coefficients as new coefficients, and the F l, s as 
additional terms. If these equations have solutions, then the solutions must be the 
partial derivatives of the F % functions. 

Now, under the assumptions of Theorem IA.21 by Theorem I A. 1 1 there is a (at 
least) continuous solution F l (x, £). Plugging that solution into the equations we just 
derived for the G"s and Ws, one obtains equations whose coefficients and boundary 
conditions are (at least) continuous. Hence Theorem IA.1I can be applied implying 
that the G l, s and H l, s are continuous. Thus F % {x, £) <E C 1 (T), proving the result. 

Appendix B. Technical results. 

Next we give some technical lemmas used throughout the paper. The first lemma 
follows from the fact that the control direct and inverse kernels are C 2 (T) functions. 
Lemma B.f. 

(B.l) 
(B.2) 
(B.3) 
(B.4) 
(B.5) 
(B.6) 

The next lemma is based on the fact that, since Anl{u, x) is twice differentiable with 
respect to u and x, and since we have Ajvl(0, x) — 0, it follows that there exists a 5a 
and K\, K2, K3 such that if |u| < 5a, then, for any v,w £ M 2 , it holds that 

dA NL (u,x) 



mi 


< 


Ci (I7I -t 


- 1 7l 


lO, 




< 


c 2 (N 4 


- 1 7l 


lO, 


I'd Ml 


< 


c 3 (l-rl 4 


- 1 7l 


l0, 


1^2 Ml 


< 




- 1 7l 


v), 


IAMI 


< 


C 5 (|7l 4 


- 1 7l 


l0, 


AiMI 


< 


Ce (|7l 4 


- 1 7l 


l0- 



\A NL (u,x)\ + 
dA NL 



da 



-(u, x)v 



dx 

dA NL (u,x) 



dudx 



o 2 a 



NL 



du 2 



[u, x)vw 



<K x \u\, 



< K 9 \v\ 



< Kz\v\\w\ 



(B.7) 
(B.8) 
(B.9) 



Similarly, since /atl(m, a;) is twice differentiable with respect to u and once with respect 
to x, and Jnl(0,x) — g L (0, x) — 0, there exists a 5f and K4, K$, Kq such that if 



\u\ < Sf, then for any v £ K 2 , 

\f NL {u,x) \ + 



NL 



dx 

OfNL 



du 



(u,x) 
(u,x) 



d 2 l 



NL 



du 2 



(u, x)v 



If 

<K B \u\, 
<K A \v\. 



Then, the following lemma holds. 

Lemma B.2. For \\j\\oo < min{<5A, 5f}, 

|Fi|<C 5 (|7l + ll7lUi), 
|F 2 |<C 6 (| 7 | 2 4-||7ll! 1 ), 
|^ 3 |<CV (1Mb || +|7l) (117,1 



L- 



|F4|<C 8 (| 7 | 2 + || 7 ||| 2 ) 



(B.10) 
(B.H) 
(B.12) 



(B.13) 
(B.14) 
(B.15) 
(B.16) 
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The next lemma follows from the previous ones. 
Lemma B.3. For ||7||oo < miri{S A ,Sf}, 



IFnl^OKM + IMUO, (B.i7) 

|Fia|<Cio(|7x| + |7l + ll7lUO, ( B -!8) 

\F 2l \ <C u (|7H-||7|Ui)(N + |NUi), (B.19) 

\F b \ < Ci 2 (|r?| + IfoM (| 7 | + ||7lU») + C14 (h| + (| 7 x| + ||7xlU») 

+C 15 |7(0)||J7(0)|, (B.20) 

|F 6 | <C 16 (M + N| i2 )(|7| + ||7lM- (B.21) 



The next lemma follows immediately from the previous lemmas and the corre- 
sponding definitions. 

Lemma B.4. For \\j\\oo < min{<$A, 8f}, 

\f 13 \ < Cvt (hi 2 + IMIlO + c ls (\e\ + \\e\\ L1 ) , (B.22) 
\f u \ < c w (hi + \\ V \\ L1 ) (1 + | 7 x| + hi + hlUO + c 20 (hxl + hi + IhlUO , (B.23) 

\F 22 \ < C 21 (| 7 | + || 7 |Ui) (\e\ + \\0\\ L i) + C 22 (h| 2 + IhlHO , (B.24) 
\F 7 \ < c 23 (hi 2 + |h||| 2 ) (1 + H7IU + ||7,||oo) + c 2i (hi + u\ L .) (h x | + H\ L ,) 
+c 25 (| 7 | + IItIU- + HtxIIoo) + \\b\W) + c 2e (h(o)| 2 + | 7 (o)||0(o)|) ,(b.25) 

|F 8 | < C 27 (h| 2 + |M|| 2 ) (1 + ||7|U) + C 28 (| 7 | + hlM + \\e\\ L 2) . (B.26) 



Finally, the following result is crucial in establishing Theorem 14. II 
PROPOSITION B.5. There exists S such that, if ||t'|| 00 + h||oo < 5, then the 
following inequalities hold 



Halloo < Cl (|| Woo + Ihlloo) j (B.27) 

\\8\\l* < C2 (hxxWv + Wlxh* + h\\v) , (B.28) 

||7xx||oo < c 3 (\\9\\oo + Ihlloo + hlloo) , (B.29) 

||7xx|U= < c 4 (\\6\\l* + \\v\\f + > (B.30) 



where Ci, 02,03,04 are positive constants. 

The proposition is proven using a series three lemmas. 

The first lemma gives a relation between the L 2 and infinity norms of 77 and j x , 
under the assumption that H7H00 is small enough. 

Lemma B.6. There exists 5 2 such that, if H7H00 < S 2 , then the following inequal- 
ities hold 



Ihlloo < Ci (|| + hlloo), (B.31) 

\\v\\l> <02(||7*IU» + ll7lUO» (B- 32 ) 

Moo <c 3 (IhHoo + hlloo), (B.33) 

||7x|U3<C4(|h|U 3 + ||7|M, (B.34) 



where 01,02,03,04 are positive constants. 
Proof. First, from (|5. 15[) we see that 



f] - S(x) 7x + F 3 [ 7 , 7 J(x) + F 4 [7](a:) = 0. 



(B.35) 
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Therefore, calling S\ the value of 5 in Lemma lB.21 and assuming H7II00 < f5i , we can 
compute a bound on h||oo as follows: 

IMIoc < ^ihxIU + Hi^^ 7»] Hoc + II-F4WII00 

< K 2 (H^lloc + II7.II00II7II00 + hllL) < K 3 (||7xlU + hlloo) • (B.36) 

Proceeding similarly with the L 2 norm, 

IMU= < ^i||7x|U= + ||-Fs[7,7*]|U» + II^HHi" 

< K 2 {\\ lx \\ L , + ||7,||i»|l7lloo + hllochlUO < K 3 (\\ lx \\ L * + 11711^X3.37) 

For the last two inequalities, we solve for j x : 

^ x = ^ 1 (x)( V + F 3 [ 7 , lx }(x)+F 4 [ 1 }(x)). (B.38) 

Remembering the definition of £(x), e = max xe [o,i] | , e3 ^ x ) } > Oj an< ^ assuming 
that H7H00 < ^1 we obtain 

||7x||oo <e (h||oo +^i||7xlU||7lloo + ^2||7llL)- (B.39) 

Therefore if we choose ||7||oo < min ^61, gi^g }; we reach the third inequality. Pro- 
ceeding similarly with the L 2 norm: 

ll7»IU» < HWvWv + tfsh.MMloo + K4<y\\v hlU) , (B.40) 

so choosing H7II00 < niin |<5i, 2 g 3 g j 1 we reac h the fourth inequality. Therefore, choos- 
ing S = min|<5i, 2 k 1 i 1 2K 3 e \> au inequalities are verified and the lemma is proven. 
□ 

The next lemma gives the relation between r] x and j xx , both in the infinity norm 
and the L 2 norm, for small ||7|| OO- 
LEMMA B.7. There exists 8 such that, *f 1 1 T 1 1 00 < S, then the following inequalities 

hold 



||7xx|U < Cl (||»fe||oo + hlloo + ||7l|co) , (B.41) 

hxxWi? < C2 {\\vx\\v + WvWi? + hlUO , (B.42) 

— C3 ( 1 1 ^fxx 1 1 00 + hlloo + hlloo), (B.43) 

WvxWl* < c 4 (hxxh* + + IItIUO - ( B -44) 

where 01,02,03,04 are positive constants. 
Proof. Taking an ^-derivative in (|5.15l) : 

j] x - E'(*)7x " Efrh** + Fi[ 7 ] 7a;a; + F 32 [j, 7 X ] + ^4 2 [7,7x] = 0, (B.45) 
where F32 and F42 are defined as: 

F 32 = JC 2 [Fx [ 7 ]7x] + F12 [7, 7x]7«, (B.46) 
F 42 = IC2 [FtlHd [7] + F 2 [j\] + Fi2[7.7*]A H 

+.F1MA1 [7] + Fi[7]i(^,a;)7x + i 71 2 3 [7,7x]7., (B.47) 
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where 

F 23 h,j x } = ^ {LbU) + ^ (£[<y],x) { lx + Ci[j\) . (B.48) 

These functionals verify the following bound, similar to the bounds developed in 
Lemma IB. 21 if ||7||co < min{(5A, Sf}: 

\F 32 \ < Cr (|| 7 || L2 + | 7 |) (||7*IU« + |7x|) + C 2 \ lx \ 2 , (B.49) 
\F 42 \ < C 4 (h\\ L 2 + | 7 |) (||7«IU« + |7x|) ■ (B.50) 

Therefore, using LcmiricL IB.6| etiid inequality (|B.13[) . and making jj'ylloo 

small enough, 

we can compute the bounds as in the proof of Lemma IB. 61 □ 

Finally, the next lemma relates rj x and 9, both in the infinity norm and the L 2 
norm, for small H7II00 and H^Hoo- 

Lemma B.8. There exists 5 such that, if H7II00 + IMIoo < then the following 
inequalities hold 

\\e\\oo < ci (Ht/xIIoo + Halloo + 1I7IU) , (B.51) 

\\6\\l*<02(\\ti x \\v + \\ti\\i* + \\'y\\l>), (B.52) 

Halloo < C 3 (H^Hoo + Nloo + IMIoo) , (B.53) 

\\rhh' < c 4 {\\0\\l* + \\vh» + II7IU0 , (B.54) 

where ci, 02,03,04 are positive constants. 

Proof. We can write (|5.29p analogously to (|5.15p : 

r] t - E(x)r) x + F 31 [7, 7x,V,Vx\{x) +F 6 [ 7 ,T}] (x) = 0, (B.55) 

where F31 is defined as: 

F 31 = K [F t [ 7 ]7fc + F n fob*] . (B.56) 

The functional F31 verifies the following bound, similar to the bounds developed in 
Lemma |B. 21 if H7II00 < rnin{5,v,5/}: 

\F 31 \ < Ci (\h\\ L 2 + | 7 |) (H^IU, + |%|) + C 2 (\\ lx \\ L 2 + | 7x |) (||77|U= + |r?|) . (B.57) 

Therefore, using Lemma IB.3[ Lemma IB.6[ and inecjuality (|B.13P , and making 1 1 7 1 1 00 ~t~ 
1 1 77 1 |oo small enough, we can compute the bounds as in the proof of Lemma lB.6l □ 
Combining the three lemmas, the proposition immediately follows. 
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